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In this paper it is proposed that the charged charmonium-like state Z c (3900) is a resonance above the threshold 
from the DTP interaction. The DTP interaction is described by the one-boson exchange model with light meson 
exchanges plus a short-range 7 /if exchange. The scattering amplitude is calculated within a Bethe-Salpeter 
equation approach and the poles near the DD * threshold are searched. In the isoscalar sector, two poles found 
under the DD * threshold, i.e., bound states, have the quantum numbers I G (J PC ) = 0 - (l + ~) and 0 + (l ++ ). The 
latter can be related to the X(3872). In the isovector sector, a bound state with I G (J P ) = 1 + (1 + ) is found with 
a large cutoff at about 3 GeV. If a cutoff at about 2 GeV is adopted with which a pole carrying the quantum 
number of the X(3872) is produced at an energy of about 3871 MeV, the pole for the bound state with 1 + (1 + ) 
runs across the threshold to a second Rienman sheet and becomes a resonance above the DD* threshold, which 
can be identified with the ZT3900). With such a cutoff, the DD * invariant mass spectrum is also investigated 
and the experimental results found by BESIII can be reproduced. 

PACS numbers: 14.40.Rt,13.75.Lb,11.10.St 


I. INTRODUCTION 

In the past decade, an amount of the so-called XYZ par¬ 
ticles were found in the facilities around the world, such as 
Belle, BABAR and BESIII. It is difficult to put these par¬ 
ticles into the conventional quark model frame, which at¬ 
tracts physicists’ special attention. An interesting observa¬ 
tion is that many of the XYZ particles are near the thresh¬ 
old of two charmed mesons. For example, the mass gap be¬ 
tween the X(3872) observed by the Belle Collaboration 0]] 
and the D°D*° threshold is smaller than 1 MeV. The struc¬ 
ture Z c (3900) observed by BESIII jD] is also very close to the 
DD* threshold. It inspired the idea that these particles origi¬ 
nate from the interaction of two mesons, such as the hadronic 
molecular state and threshold effect. 

In the literature, many efforts have been made to study the 
possibility of interpreting the X(3872) and the Z, (3900) as 
a hadronic molecular stat (that is, a bound state below the 
threshold), such as calculations with the QCD sum rule yl- 
H. In Refs. SH, the BB*/DD* system was studied with a 
nonrelativistic one-boson exchange (OBE) model by solving 
the Schrodinger equation. There does not exist a bound state 
from the DD* interaction, which can be identified with the 
Z c ( 3900) observed by BESIII. In Ref. fst], the DD* interac¬ 
tion is studied in a Bethe-Salpeter equation approach. A state 
with quantum number I G (J P ) = 0 + (l ++ ) is produced from the 
DD* interaction, which corresponds to the isoscalar particle 
X(33,12). No bound state related to the Z c (3900) was found. 
In lattice calculations, a candidate X(3872) state was observed 
while the possibility of a shallow bound state related to the 
Z c ( 3900) was not supported JM1- 

Generally, the theoretical studies suggest that there exists a 
bound state relevant to the AT3872) from the DD* interaction, 
while the existence of a bound state relevant to the Z c ( 3900) 
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is disfavored. Besides the tetraquark interpretation l fl2l - [l4ll , 
many authors proposed an alternative explanation that the 
structure Z f (3900) is simply a kinematical effect 13Hi- In 
their opinion, such structure is not related to an S-matrix pole 
and therefore should not be interpreted as a state. In Ref. 
d, it was suggested that in the elastic channel the kinematic 
threshold cusp cannot produce a narrow peak in the invari¬ 
ant mass distribution in contrast with a genuine S-matrix pole, 
which can be used to distinguish kinematic cusp effects from 
genuine poles. 

In scattering theory, a peak structure in the experiment can 
be related to not only a bound state below threshold but also 
a resonance above threshold. Both the bound state and reso¬ 
nance are from an attractive interaction, but the former needs 
stronger attraction 0. For example, a popular interpreta¬ 
tion of the A(1405) is a dynamically generated state with a 
two-pole structure lfl9U2lll . It is interesting to note that the 
higher-energy KN channel has a stronger attraction to support 
a bound state, while the lower energy ttY. channel shows a rel¬ 
atively weaker attraction, which is nevertheless strong enough 
to generate a resonance 12211 . In all measured channels, the ex¬ 
perimental mass of the Z c (3900) is higher than the threshold 
1.2311, which also indicates that it may be a resonance instead 
of a bound state. Hence, it is interesting to study the possibil¬ 
ity of interpreting the Z r f3900) as a resonance above the DD* 
threshold instead of a bound state. 

In the literature, study about the possibilities of the XYZ 
particles as a resonance has been scarce. There exist many 
works to study the bound state from meson-meson interac¬ 
tions by the potential model, the QCD sum rule and the lat¬ 
tice calculation. However, it is difficult to deal with a reso¬ 
nance with available QCD sum rule or lattice technology. In 
Ref. J3, the poles of a T matrix were searched in the lo¬ 
cal hidden gauge approach with heavy quark spin symmetry, 
but only the bound state from the DD* interaction was stud¬ 
ied. In this work, the method in Ref. il will be extended to 
search the poles for both bound states and resonance, and the 
heavy quark effective theory will be used to describe the DD* 



2 


interaction with light meson exchanges plus short-range J/ip 
exchange. 

This work is organized as follows: In the next section a 
theoretical frame is developed based on a quasipotential ap¬ 
proximation of the Bethe-Salpeter equation. In Sec. Ill, the 
potential kernel with light meson exchange and J/if/ exchange 
is derived with the help of the effective Lagrangian from the 
heavy quark effective theory. The numerical results are given 
in Sec IV. In the last section, a summary is given. 

II. SCATTERING AMPLITUDE 

The scattering amplitude can be obtained through solving 
the Bethe-Salpeter equation. The general form of the Bethe- 
Salpeter equation for the scattering amplitude as shown in 
Fig. Q] reads, 

M(k\k’ 2 ,k i k2-,P) 

/ d 4 k" 

■ 'V{k\k' 2 ,k”k 2 \ P)G(k';k")M(k"k",k\k 2 \P), (1) 

where 'V is the potential kernel and G is the product of the 
propagators for two constituent particles. Here the momentum 
of the system P - k\+k 2 = k\ +k 2 , and the relative momentum 
k" = (k" - k")/2. 



FIG. 1: The diagram for the Bethe-Salpeter equaiton. The thick and 
thin lines are for particles 1 and 2, respectively. 


Note that the sum extends only over non-negative A", and a 
factor 1 / V2 has been included in the scattering amplitude and 
potential for zero helicity. The potential is defined as 

'V&fe'.p) = 2tt Jd cos 6 [d J u ,(0yV A , A (p\p) 

+ (4) 

where the momenta are chosen as k\ = (W - E, 0,0, -p), 
k 2 = (£, 0,0, p) and k' x — (W - E', —p' sin 6, 0, —p' cos 6), 
k' 2 = (E',p' sin$, 0,p' cos9) with p = \p\ in order to avoid 
confusion with the four-momentum p. 

The above equation can be related to the Lippmann- 
Schwinger equations used by Oset et al., if the potential ker¬ 
nel f V is only dependent on the square of the momentum of 
the system s = P 2 and the Go is chosen as the one used in 
Ref. |[2Ct l. A cutoff regularization has been introduced in the 
integration of the propagator in Ref. (20!], and it is related to 
a dimensional regularization. Since in our formalism, the in¬ 
tegration is on the potential also the only cutoff regularization 
is practical here. In this work we will adopt an exponential 
regularization instead of cutoff regularization by introducing 
a form factor in the propagator as 

G 0 (p) -> G 0 (p) [ e - ( *?-m]) 2/A4 1 2 . ( 5 ) 

Here particle 2 is not involved in the form factor due to its 
on-shell-ness. The exponential regularization used here can 
be seen as a softer version of the cutoff regularization in Ref. 
0, where the momentum p is cutoff at a certain value p max- 
The cutoff A plays an analogous role to the cutoff p max of cut¬ 
off regularization. It can also be understood as a form factor 
in exponential form for the charmed mesons to reflect the in¬ 
ternal structure of the hadron and to make the integration con¬ 
vergent. It is consistent with the OBE model where a form 
factor is usually added for the off-shell particle. 


In this work, a quasipotential approximation will be applied 
to reduce the four-dimensional Bethe-Salpeter equation to a 
three-dimensional one. Here the covariant spectator theory 
will be applied as shown in Appendix A, in which the heav¬ 
ier meson, particle 2, is put on shell. After multiplying the 
polarized vector on both sides of the equation, we have 


III. THE POTENTIAL 

In Ref. (H], it was explained explicitly how to construct a 
potential for states with definite isospin under SU(3) symme¬ 
try with the corresponding flavor wave functions [g] 


Ma’Ap\p) = *v xx (P'-P) + Zf yy 


'Vw( P ',p")G 0 { P ")Ma"Ap”,p), ( 2 ) 


where p, p' and p" are the momenta of constituent 2. And the 
potential e V xx ip', p) = e l "(p'fV llv (p',p)e]' l (p). 

To reduce the Bethe-Salpeter equation to a one-dimensional 
equation, we apply the partial wave expansion as shown in 
Appendix B. The partial wave Bethe-Salpeter equation with 
fixed parity reads as, 


Kxi p'.p) - tUp'.p) + 


y r P"W 

V J W 


Vi, A „(p',p")Go(p")Ml„ A (p”,p). (3) 


|Z+y /=1 = -^=(\D* + D°) + c\D + D*°)), 

I z Db ,)i=\ = -^(|ZJ*-5 0 ) + c|Z)-5* 0 », 

= \[(\D* + D~)-\D t0 D°)) 

+ c(\D + D*-) - |£>°5*°»], 

I z "p,)/=o - i[(|o* + zr> + |D*°fi°» 

+ c(|Z) + £)*~) + |Z)°Z)* 0 »], (6) 

where c = ± corresponds to C-parity C - + respectively. For 
the isovector state the c is related to the G-parity. 









3 


Basically, the strong interaction should be described by 
gluon and quark freedoms. If the distance between two 
hadrons is large, the interaction will appear as meson ex¬ 
change, that is, the OBE model. Many efforts have been made 
to study the connection between QCD and the OBE model. 
For example, in Refs. BEi after integrating out quark 
degrees of freedom in the effective Lagrangian, one obtains 
the chiral effective Lagrangian for mesons. Though there is 
still not particularly convincing about the connection, as a 
phenomenological model the people’s confidence about the 
OBE model arises from the successes of its_ applications to 
the deuteron, such as the CD-Bonn model [27] and the Gross 
model 128 [, and the constituent quark model by Riska and 
Glozman 1291. If the X(3872) is interpreted as a molecular 
state, it should have a radius of about 7 fm, as estimated by 
Close and Page lEfJ. For a resonance above the threshold, it 
is reasonable to assume the interaction is at a large distance. 
Hence, the OBE model will be adopted to describe the DD* 
interaction in this work. 

It seems strange to include the vector-meson exchanges, 
which mainly take effect at short distances. However, we 
should remind the reader that a short-distance interaction does 
not mean an interaction only at short distances. The vector- 
meson exchange will have some remnant at long distances. 
It is also can be understood as indicating that there is some 
possibility that the two hadrons are close to each other. If the 
interaction from the vector-meson exchange is large enough, 
it will have a considerable remnant at long distances. Based 
on this consideration, the vector-meson exchanges are also in¬ 
cluded in the OBE model. At short distances, the studies in the 
constituent quark model suggests the gluon exchange contri¬ 
bution can be replaced by the contribution from vector-meson 
exchanges 113 lL l32ll . Therefore, vector-meson exchanges will 
be included in the calculation in this work instead of gluon 
exchange. 

There exist two types of diagram namely, the direct dia¬ 
gram and the cross diagram (see Fig. [2}, as in a conventional 
OBE potential model 0, ..71 ■ In the cross diagram, the final 
particles are alternated. With such alternation, the propagator 
is the same for different components in a SU(3) state, such as 
| D* + D {) ) and \D + D*°) for the positive-charge state, so that the 
equations for different components are reduced to one equa¬ 
tion under S(3) symmetry [|8|]. 


k 


I 
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FIG. 2: The potential including both direct and cross diagrams. I d 
and 1^ are the isospin factors for direct and cross diagrams, respec¬ 
tively. 


To write the potential, we adopt the effective Lagrangians of 
the pseudoscalar and vector mesons with h eavy flavor mesons 
from the heavy quark effective theory 1133113411 . The potential 
kernels *V from vector-meson (V) exchange, pseudoscalar- 
meson (P) exchange, and scalar-meson (cr) exchange have 


been given in Ref. 0 as 


r\/Direct 

K V 


r 


•Direct 


<v- 


Cross 
V 


rx/Cross 


,P 2 g 2 v {ki + K ) • {k\ + k\) ei ■ e' 2 

2 q 2 ~ 

.. 2 62 ' e 2 

i4g lr m r mr- —--y, 

q A — m£. 

i2A 2 g 2 v £A a pii(k2 + k\Yq a <% 

EA'a'pvikx + k^Y'q a ' —r —— 
q- - m * 

. 4 g 2 mpmp* k’ x ■ 62^1 • 
f 2 q 2 - ml 


(7) 


where the momenta k\ \ and q are defined as in Fig. [2] Here 

The 4 1 is the polarization vector for the initial or final parti¬ 
cle 2. m p, my and m ir are the masses for the exchange pseu¬ 
doscalar, vector, and sigma mesons. In the OBE model, the 
masses of all mesons used are from the PDG [[23]. The cr me¬ 
son is seen as a particle with a mass 475 MeV, which is the 
center value suggested by the PDG. 

In the OBE model, we will adopt the physical values of the 
coupling constants. How to determine the coupling constants 
has been discussed in the heavy quark effective theory & 
The coupling constant g for the pseudoscalar exchange is ex¬ 
tracted from the experimental width of D* + as g - 0.59 f35i l. 
The parameter fi for the vector meson was fixed as (A = 0.9 by 
the vector-meson dominance mechanism and A = 0.56 GeV 1 
was obtained by comparing the form factor obtained by lattice 
QCD with the one calculated by the light-cone sum rule [330- 
The coupling constant go- - g „/(2 V6) with g n — 3.73 was 
given in Ref. iB l. 

A form factor is introduced to compensate the off-shell ef¬ 
fect of the exchange meson as f(q 2 ) = ( A 4 g2 ) 4 - It is different 
from the one used in the propagator [see Eq. ©] , which is 
usually used in a case where the off-shell particle has pos¬ 
sibility q 2 = A 2 , to avoid an unnecessary pole arising from 
the form factor. The cutoff can be related to the radius of the 
hadron as 


r 


2 


6 d 
F( 0) dq 2 


F(q 2 )\ q 2 =0 . 


( 8 ) 


With Eq. ©I, the cutoffs in two types of form factors have 
relation A « 5/r. If we assume the mesons have radii of about 
0.5 fm, the cutoff is about 2 GeV. Here the momenta for the 
exchange mesons are defined as q - k' 2 - ki and q — fc' - ki 
for direct and cross diagrams, respectively. In the propagator 
of the meson exchange we make a replacement q 2 —> —| q 2 \ to 
remove the unphysical singularities, as in Ref. 113711 . 

In Ref. I 24ll . it was suggested that the / /if/ exchange is im¬ 
portant in the DD* interaction. And the potential was writ¬ 
ten in the local hidden gauge approach with heavy quark spin 
symmetry. In this work, the couplings of heavy-light charmed 
mesons to Jlip are written with help of the heavy quark effec- 
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tive theory as J34ll38l] 


£»WI# = ' O* d D * 

- i (SD*-*d p D*+ i// p D*-*d D* p )\, 

-Cd U d< s) ji* = g Dls) D^ e^ aT cfriDd T D m + D* a ^d T D), 

_ ^^ 

-CduD^jw = igD (s) D^iff D d D, (9) 

where the couplings are related to a single parameter g 2 as 

gD*D*i// gD (s) D^ - 

-=-= ,?»• ,«(,# = 2g 2 V^, (10) 

»l£> W ' 

with g 2 = y/mj;/(2m D fy) with = 405 MeV. 

With the above Lagrangians, the potential kernel for J/ijj 
exchange is written as, 

^/r = -igD-D'JItl/gDDJIfi^2 ' (M + &i) €2 ' (^2 + 

+ e'2-{k 2 + k' 2 ) e 2 -{k x +k\) 

- ( k 2 + k' 2 ) ■ {k\ + k\) e' 2 ■ e2 ]—- - —, 

r - m m 

^/T = ig 2 DD-Ji^ E ^(k2+k\Yq a ^ 

■ SA'a'/3v(k\ + k' 2 y V q a 'e '2 — — 2 —■ (11) 

<r - m J!t 


We would like to point out that the potentials obtained by 
the heavy quark effective theory are comparable to the ones 
obtained from the chiral Lagrangian in Ref. [24] after a non- 
relativization. For example, the potential kernel for the J/ijj 
exchange can be rewritten as 


rx/Direct 

v jo* 


gD'D‘J/tpgDDJ/1/ 

2 

m T/ . 

-W 

„2 _ • 


(*2 + k' 2 ) -{h+k[)e' 2 - 62 ,( 12 ) 


with g D >D‘JltgDDJii,lm 1 Jliji = 7.04 in this work and 


TABLE I: The isospin factors If and If for direct and cross diagrams 
and different exchange mesons. 



Direct diagram 

Crossed diagram 


p U) J / l/r cr 

71 Tj p Cl) 

J/> A 

[prf 

---11 

2 2 x 1 

-i c i c - l -c l -c 

2 c 2 c 2° 

C 

\pp , ] s 

1 i 1 1 

2 2 

lr i c lr ±C 

2 c 6 c 2 C 2 c 

C 


A. Numerical solution of the Bethe-Salpeter equation 


Before solving the one-dimensional partial wave Bethe- 
Salpeter equation numerically, we need to deal with the pole 
in Go(p) as (here the notation J p is omitted) 


/M(p.p') = rV(p, p') 

p" 2 dp ". 


f 


(2;r) 3 


i‘V(p,p")Go(p")iM(p".p') 
dp" A(p") . ,p" 2 <5Go(p") 


mP ’ P::)[ J (2?r) 3 p" 2 -p"2 

iM(p",p')9(s -nu — m 2 ), 


+ l “ 


8n 2 


(14) 


with 


n" 2 

A(Vo) = [p'V ,2 -p" 2 )Go(p")W' = 


tfG 0 (p , W'-p") = «y(G„(p) ,, ) = 


1 


4Wp" 


d(p"-p"), (15) 


with p" = ^ yj[W 2 - (M + m) 2 ][W 2 - (M - m ) 2 ]. 

To solve the integral equation, we discrete the momenta p, 
p', and p" by the Gauss quadrature with wight w{ p,) and have 


<v, j; r ; ct = ~(k 2 +k’ 2 )-(k [ +k\) £ ’ 2 -e 2 , ( 13 ) 

with Cij/Af 2 = 6.89 in Ref. iHl . The main difference be¬ 
tween our work and Ref. |24ll is that a form factor is added not 
only to the light pseudoscalar-meson exchanges but also to the 
J/ifi exchange, which will suppress the contribution from J/ip 
exchange. 

The flavor factors if and If for direct and cross diagrams are 
presented in Table Q] The cancellation of the p and at meson 
exchanges happens in the isovector sector as suggested by the 
isospin factor listed in Table [T] which leads to a shortage of 
the short-range interaction in that sector if the J/fj exchange 
is absent. 


IV. NUMERICAL RESULTS 

In this work, we will search the poles of the scattering am¬ 
plitude A1 from the DD* interaction which is described by the 
potential kernel obtained in the above section. The scattering 
amplitude is obtained by solving the Bethe-Salpeter equation. 


N 

M ik - V ik + J]VijGjMj k , (16) 

2=0 


where i is absorbed in M or V. The discreted propagator is of 
a form 


G/>o 
G i =0 


w( P p P f 

(2tt) 3 


Co(p"), 


‘Po + y MPj) P" 2 
327 t 2 W y (27 t) 3 2 W(p" 2 - p" 2 ) 


(17) 


In this work, we will search the poles from the amplitude 
of elastic scattering where the initial and final particles are on 
shell. The scattering amplitude is 

M = Moo = J][(l - VGf-'hjVjo. (18) 

i 


The pole can be searched by variation of 2 to satisfy 


|l-V(z)G(z)[ = 0, (19) 
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where z - Er + iT/2 equals the meson-baryon energy W at 
the real axis. Since z — + p 2 + „ Jm\ + p 2 , the p-plane 

corresponds to two Reimann sheets for z. A bound state is 
located in the first Reimann sheet while a resonance is located 
in the second Reimann sheet with lm(p)<0. Since only one 
channel is considered in this work, the bound state is located 
at the real axis of the z complex plane, while the resonance 
will deflect the real axis to the complex plane. 

B. Bound state from 1 ) 1 )' interaction with a A scan 

There exist two types of pole, namely bound state and res¬ 
onance. First, we will make a A scan from 0.8 GeV to 4 GeV 
to find the bound state from the DD* interaction. The pole of 
a bound state from DD* interaction is located at the real axis. 

TABLE II: The position of the bound state from the DD * interac¬ 
tion at the real axis Re{z) = W with a A scan. The second and third 
columns are for the full model. The results without the y/i/r exchange 
are listed in the fourth and fifth columns and compared with the re¬ 
sults in Ref. fl]. The JI\Ji (I) in the eighth and ninth columns and 
7/ 1 A (II) in tenth and eleventh columns are for the results from only 
J/i// exchange with and without form factor. The cutoff A and energy 
W are in units of GeV. 



Full model 

No 7/l/r 

Ref. [8J 

W (I) 

J I* (II) 

I G (J PC ) 

A 

W 

A 

W 

A 

W 

A 


A 

W 

0-(0“) 











+ 

1 

O 

+ 

O 











o-(i— ) 











0 + (l- + ) 











o-(V-) 

1.0 

3.864 

1.0 

3.868 

1.3 

3.876 

- 

- 

2.5 

3.867 


1.2 

3.848 

1.2 

3.854 

1.4 

3.870 

- 

- 

2.6 

3.850 

0 + (l ++ ) 

1.9 

3.873 

1.9 

3.875 

2.0 

3.876 

- 

- 

3.9 

3.875 


2.4 

3.871 

2.4 

3.874 

2.4 

3.872 

- 

- 

4.0 

3.836 

U(O-) 











V(O-) 











U(l-) 











1(1) 











U(1 + ) 

3.0 

3.874 

- 

- 

- 

- 

- 

- 

2.4 

3.875 


3.3 

3.858 

- 

- 

- 

- 

- 

- 

2.5 

3.867 

1-(1 + ) 












One can find that the //i/r exchange plays a more impor¬ 
tant role in the isoscalar sector than in the isovector sector, 
where the short-range interaction is absent due to the cancel¬ 
lation between p and at exchanges in the isoscalar sector, as 


shown in Table |T] As expected, the results without the ex¬ 
change are also close to those obtained from the solution of a 
Bethe-Salpeter equation for the vertex in Ref. [01 . The results 
with only the J/i// exchange are also listed in Table HU There 
is no bound state found for the J/ip exchange with form factor 
[labeled as 7/<A (I)], while if the form factor is removed [la¬ 
beled as J/ifj (II)], the bound state is found as in the full model, 
which indicates the form factor weakens the contribution from 
the J/ip exchange. 

In the isoscalar vector, there exist hidden charmed bound 
states with I C (J PC ) = CL(1 + ~) and CL(1 ++ ) with cutoffs about 
1 GeV and 2 GeV, respectively. The DD* bound state with 
0 + (l ++ ) can be related to the A(3872). In the isovector sector, 
a bound state with 1 + (1 + ) can be found with a larger cutoff at 
about 3 GeV. 


C. The Z c (3900) as a resonance 

In physics, the cutoff in the DD* interaction should be same 
for different quantum numbers. Usually, a decrease of the 
cutoff will lead to a weaker interaction, and vice versa. As 
stated in the scattering theory, when interaction weakens, the 
bound state runs to the threshold and becomes a resonance 
if the interaction is still strong enough. If the cutoff is in¬ 
creased to 3 GeV, the pole with 0 + (l ++ ) will move to about 
3650 MeV, which is very far from the experimental mass of 
V(3872). Hence, we decrease the cutoff A for 1 + (1 + ) from 
about 3 GeV to 2.4 GeV with which a bound state which 
has quantum number 0 + (l ++ ) is found at 3.871 GeV. A pole 
with 1 + (1 + ) is produced slightly higher than the DD* thresh¬ 
old with a cutoff A = 2.4 GeV as shown in the upper panel of 
Fig. [3 It is located at z = 3876 + ;5 MeV and can be identified 
with the charged charmonium-like state Z r (3900) observed in 
BESIII. 

Obviously, the position of the pole is below the experiment 
masses ofZ c .(3900), 3883.9+1.5+4.2MeV in the DD* channel 
and 3899 ± 3.6 ± 4.9 MeV in the channel fzill . The 

experimental mass is obtained by fitting the invariant mass 
spectrum. As shown in the middle and lower panels of Fig. [3 
the invariant mass spectrum of the DD* channel is presented 
and compared with the experimental results released by the 
BESIII Collaboration HJ. 

The invariant mass distribution is given approximately as 

Hi 

— = C|M 1+ | 2 T3 (W 2 , M 2 , m 2 )A^ (W 2 , W 1 , ml)/W, (20) 
dW 

where M 1+ is the scattering amplitude obtained from the 
Bethe-Salpeter eqaution as defined in Eq. (IT8] >. W is total en¬ 
ergy of the process, and in-] is the third final particle n J24, 41]. 
A background spectrum is also included as in Refs. 12413911 . 

fbk g (W) = a(W - M min ) c (M max - W) d , (21) 

where M mm and M max are the minimum and maximum kine¬ 
matically allowed masses. The general constant C in Eq. (f20b . 
and the general constant a and exponents c and d in Eq. ([2Tb 
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FIG. 3: (color online) The ]1 - G(z)V(z)| for I G (J PC ) = 1 + (1 + ) in the 
complex energy plane (upper panel) and the D°D*~ (middle panel) 
and D + D*° (lower panel) invariant mass spectra. The abscissa axis 
Re(z ) represents the corresponding DD * invariant mass W in units of 
GeV. The solid, dotted, and dashed lines are for results with the full 
model, the resonance contribution, and the background, respectively. 
The experimental results are from BESIII [1391 .The theoretical results 
are normalized to the experiment. 


are free parameters adjusted to reproduce the experimental 
data. 

The experimental data are well reproduced with the reso¬ 
nance contribution plus a background. A peak is found at 
3881 MeV which is higher than the position of the pole of the 
resonance but closer to the experimental results 3883.9+1.5 ± 
4.2 MeV in the DD* channel l lvill . We would like to note 
that due to the pole’s nearness to the threshold the contribu¬ 
tion from the Z e (3900) in our model does not have a standard 
Breit-Wigner form, which was adopted in the experimental 
fitting |39il ■ The nearness also results in a relatively large tail 
of the resonance at higher energies. These differences lead to 
a background contribution, which is quite different from the 
one in Ref. [[39] to be used to reproduce the invariant mass 
spectrum. 


V. SUMMARY 


In this work, the DD* interaction is studied within the OBE 
model, including the contribution from light meson exchanges 
plus the short-range J/ijj exchange. The scattering amplitude 
is calculated within a Bethe-Salpeter equation approach, and 


the poles near the DD* threshold are searched. It is found 
that the charged charmonium-like state Z, (3900) can be inter¬ 
preted a resonance above the threshold from the DD* interac¬ 
tion. 

In the isoscalar sector, the poles are found under the DD* 
threshold and have the quantum numbers I G (J PC ) - 0”(l + ~) 
and 0 + (l ++ ). If the J/ip exchange is excluded, the position 
of the pole is almost unaffected. In the isovector sector, where 
the short-range contributions from u> and p exchanges are can¬ 
celed, a bound state is found with I G (J P ) = 1 + (1 + ). It disap¬ 
pears if the J/ifj exchange is removed. The results show that 
the short-range J/ip exchange, which is not included in the 
conventional one-boson exchange model, is important to pro¬ 
vide an attractive interaction to produce the pole in the isovec¬ 
tor sector. 

If a cutoff A = 2.4 GeV is adopted with which a pole car¬ 
rying the quantum number of the X(3872) is produced at an 
energy of about 3871 MeV, the pole for the bound state with 
1 + (1 + ) runs across the threshold to a second Rienman sheet 
and becomes a resonance at 3876 + ;5 MeV, which can be 
identified with the Z, (3900). The line shape of the invariant 
mass spectrum in the DD* channel is also investigated and 
the experimental results by the BESIII Collaboration can be 
reproduced. A peak is found in the DD* invariant mass spec¬ 
trum at about 3881 MeV, which is higher than the resonance 
pole but closer to the experimental values. 

Note: Just after the paper was submitted, we noticed a re¬ 
lated work released in arXiv, in which the author also found a 
second-sheet pole in studying elastic DD* scattering j!42]. 
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Appendix A: The quasipotential approximation of the 
Bethe-Salpeter equation 

It is popular to reduce the Bethe-Salpeter equation from 
a four-dimensional integral equation to a three-dimensional 
equation by quasipotential approximation. In principle, in¬ 
finite choices can be applied to make the quasipotential ap¬ 
proximation. The popular methods used in literature include 
the BSLT approximation, the K-matrix method, the instan¬ 
taneous approximation, and the covariant spectator theory 
(CST) 0123ESS- 

The reduced propagator G under quasipotential approxima¬ 
tions should satisfy the unitary condition, i.e., the relation 

G - G ] = 2m6((r l] (s)P + kf - m 2 )6((ii 2 (s)P - kf - m 2 ),( Al) 

where k\ — rj\(s)P+k, k 2 — t] 2 (s)—k, and m\ 2 are the momenta 
and mass of constituents 1 and 2 with rji(s) + rj 2 (s) - 1 and 
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s = P 2 . One can define 77^2 = ei, 2/(^1 + e 2 ) with ei^s) = 
(s + m 2 2 - m 2 j)/2 \[s. Now we have many choices to write 
the propagator. The most popular form of propagator is E§- 

m 


G - 2n 


f , — ~h(s',s) 

J s' - s + le 

2 


6{W x {s')P’ + kf - m]) 6 ([tj' 2 (s')P' - kf - m 2 ),(A2) 


where P' — s/s'/sP. It is random to some extent to choose 
h(s', s) . The h function is chosen as h(s' - s) = 1 with 
rj'(s') - rj(s') for the BSLT formalisms,. In the CST, 
h(s',s) = (^fs'+ yfs)/ sis', with//j(s') = r/\(s)y/s/s' and 
r]' 2 (s') = 1 - pi(s) sfsjs 7 . 

Eqation ( IA2b is quite far from exhausting all possible three- 
dimensional reductions. For example, the widely used instan¬ 
taneous approximation has 


G = 


-f 


-1 


= in- 


° (ki - mi + ie)(jf - m ; + ie ) 
( E\(p) + E 2 (p))/E l (p)E 2 (p) 
W 2 -(Edp) + E 2 (p)) 2 ’ 


(A3) 


which satisfies the unitary condition also. 

As in Ref. @1, the covariant spectator theory 0 S 
is adopted to make the quasipotential approximation to re¬ 
duce the four-dimensional Bethe-Salpeter equation to a three- 
dimensional equation in the current work. Written down in 
the center-of-mass frame where P = (W. 0), the propagator in 
the CST is 


G = 2ni 


= 2m 


S + (k 2 


■ m\) 




6 + (k° 2 - E 2 (p)) 


2 E 2 (p)[(W - E 2 (p)f - E 2 Ap)f 


(A4) 


where k\ - (kf -p) - (W - E 2 (p),-p) andk 2 = (k 2 ,p) = 
( E 2 (p),p ) with E\ 2 ip) = yjm 2 o + \p\ 2 - In our formalism a 
definition Go = G!2n is used for convenience. 


Appendix B: The partial wave expansion and the amplitudes 
with fixed parity 


where the momenta are chosen as k\ = (W - E, 0,0, -p), 
k 2 = (£, 0,0, p) and kj = (W - E', —p' sin 6, 0, -p' cos 6), 
k' 2 — (E/.p'sin0,0,p'cost?) with p = \p\ in order to avoid 
confusion with the four-momentum p. The particle helicities 
A are the projections of the spin s on the direction of motion 
of the particle. Here and hereafter, the individual helicities are 
omitted where redundant and the states are only labeled by the 
total helicities d. A' and A". Thus, once in the center-of-mass 
system the z-axis is chosen along the three-momentum of the 
incoming particle 1, and one has A\ - Si for final state particle 
1 and A 2 = -s 2 for final state particle 2. 

The amplitudes with definite parity can be constructed as 

JH 


= M'n+TiMi^ (B4) 

where // = PP\P 2 (-\) J ~ ] 'with P and P\ 2 being the pari¬ 
ties and J and J 1 .2 being the angular momenta for the system 
and particle 1 or 2. It is easy to find that the amplitudes with 
definite parity have properties such as 

<2 = M J Zu = TfM*. (B5) 

The potential r V J A , A has analogous relations. 

The Bethe-salpeter equation for definite parity can be writ¬ 
ten as 






(B6) 


Please note that there exists a factor 4 on the second term in 
the right side of the equation. 

By using the relation in Eq. (IB5b . Eq. (lB6b can be rewritten 
as 

m j a , a = v j ; A , + ^ 'v j ; a „gm j ; u , ( B7 ) 

A">0 

where A, A' and A" > 0 and M J A , X = fA'f\M J A , v with /o = -^= 
and f A+ o = 1. 

The sum of the square of the amplitude can be written as a 
from with definite parity as: 


The partial wave expansion of the scattering amplitude M 
in Eq. (J2]» is 01 

M A ’ A (p',p) = Ji 0', °rEt'2,2 ff (P'> P) 

JA r 

■ D ] Ar A ^,6, 0), (Bl) 

where J is the angular momentum for the partial wave consid¬ 
ered and />j* ( (7/>, 6 ,0) is the rotation matrix with Ar being the 
helicity of the bound state. The potential in the partial wave 
expansion equation (0 is 

tw./o - Z P ) 

JAr 

■ D J lR ^,e, 0). (B2) 

Without loss of the generality, we choose the scattering to be 
in the xz plane, the potential is written as 

< V J iU (p,p') = 2n J d cos Qd J AA ,(G P ’ tP )'V A ’ A (p' ,p), (B3) 


Zl4/= Z !^I 2 - ( B8 > 

J,A'A J p ,A'>0A>0 

By using the normalization of the Wigner D matrix, the in¬ 
tegration of the amplitude is 

dD.\MA’A(p\p)\ 2 = Z l^(P'.P)l 2 - (B9) 

J P ,A'A 

Since there is no interference between the contributions from 
different partial waves, the total cross section can also be di¬ 
vided into partial-wave cross sections. Since only the square 
of the amplitude is related to the physical observables, we 
omit all hat notation (") if not necessary and keep in mind 
that there is a factor 1 / a/ 2 in the potential also. 
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